Morphogenesis of growing soft tissues 
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Recently, much attention has been given to a noteworthy property of some soft tissues: their 
ability to grow. Many attempts have been made to model this behaviour in biology, chemistry and 
physics. Using the theory of finite elasticity, Rodriguez has postulated a multiplicative decomposi- 
tion of the geometric deformation gradient into a growth-induced part and an elastic one needed to 
ensure compatibility of the body. In order to fully explore the consequences of this hypothesis, the 
equations describing thin elastic objects under finite growth are derived. Under appropriate scaling 
assumptions for the growth rates, the proposed model is of the Foppl-von Karman type. As an 
illustration, the circumferential growth of a free hyperelastic disk is studied. 

PACS numbers: 87.17.Rr, 46.32.+X, 46.70.Hg, 87.17.Pq 



Biological tissues are conventionally classified into two 
categories: hard tissues (e.g. bones or teeth) and soft 
tissues (e.g. muscles, arteries, tendons, skin), depending 
on their mechanical properties. Soft tissues, which typ- 
ically exhibit anisotropic, nonlinear, inhomogeneous be- 
haviours, are often subject to large stresses and strains. 
The theory of finite elasticity therefore forms an appro- 
priate framework to describe their properties [l|, 0, , 
in the absence of viscolelastic efi'ects. Along these lines, 
much work has been done to establish constitutive re- 
lationships for specific biological materials such as the 
skin, blood vessels, lung, brain, liver and kidney 0, 
although computing stresses and strains under applied 
external loads remains a difficult task. 

Observation of biological tissues has revealed the exis- 
tence of internal stresses, even in the absence of external 
loads. These residual stresses are induced by growth Q 
and affect the geometrical properties of tissues. Soft tis- 
sues may undergo volumetric growth, that is a smooth 
volumetric resorption of the bulk material [1, Q depend- 
ing on space, orientation and the state of stress within 
the body. Growth is a complex process involving bio- 
chemical and physical reactions at many different length- 
and time- scales, that occur through cell division, cell en- 
largement, secretion of extracellular matrix or accretion 
at surfaces. The removal of mass is referred to as atrophy 
and occurs through cell death, cell shrinkage or resorp- 
tion. The general idea is that the total deformation of 
the body can be due to both change of mass and elastic 
deformations 0, i, S, E, [ill, [H ■ 



Before (resp. after) the deformation, the body is in 
the reference (resp. current) configuration and the place 
of each material point is denoted by X (resp. x). We 

define the geometric deformation tensor by F = to 

describe locally the overall deformation process. In order 
to model the growth process, we follow Rodriguez et al 
[l3| in making the following three assumptions: (i) there 
exists a zero-stress reference state; (ii) the geometric de- 



formation gradient F admits a multiplicative decompo- 
sition of the form F = AG where G is a growth tensor 
describing the change in mass and A an elastic tensor 
characterizing the reorganization of the body needed to 
ensure compatibility (no overlap) and integrity (no cav- 
itation) of the body; (Hi) the response function of the 
material depends only on the elastic part of the total 
deformation. Despite its simplicity, Rodriguez theory is 
yet to be investigated, because of the complexity of fi- 
nite elasticity although inhomogeneous and anisotropic 
growth has been studied in details in some simple geom- 
etry , • In order to explore the consequences of Ro- 
driguez's model, we have written the equilibrium equa- 
tions of a thin elastic body subject to growth-induced 
finite displacements. This reduction of dimensionality al- 
lows deep theoretical inspection and is relevant for many 
biological systems (e.g. leaves, skin). 




FIG. 1: A few examples of buckling in nature. (A) A potato chip 
adopts a saddle shape during frying. (B) Acetabularia schenckii: 
a green algae. (C) A dead leaf. (D) A leaf infected by the Cotton 
Leaf Crumple Virus (CLCrV), known to affect growth and induce 
curling and the appearance of blisters at the surface of the leaf. 



Under appropriate scaling assumptions, the resulting 
equations are found to be an extension of the well known 
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Foppl von Karman (FvK) model, a powerful theory for 
buckling instabilities, that are widely diffused in nature 
(FIG. [1]), but which is also able to explain complex post- 
buckling phenomena such as crumpling. Experimentally, 
it has been shown that growth may affect curvature in 
various systems. In growing gels, both homogeneous 
growth under constraints [l6| and free inhomogencous 
growth [l3| have been investigated. Thermal expansion 
(FIG.l-A), as well as dessication, can also bend an elas- 
tic body and cause it to crumple as seen in dead leaves 
(FIG.l-C). In living tissues, viruses such as the Cotton 
Leaf Crumple Virus (CLCrV) modify the growth pro- 
cess and infected plants exhibit curled or crumpled leaves 
(FIG.l-D) but buckling can also occur during normal de- 
velopment. Some mushrooms' or algae's caps (FIG.l-B) 
may undergo symmetry breaking, and adopt an oscilla- 
tory or cup shape. At the cellular level, a new milestone 
was reached with the discovery of the CINCINNATA 
gene whose local expression affects growth and curvatures 
of the Antirrhinum (snapdragon) leaf 18| . Complemen- 
tary to the inhomogeneity of growth, anisotropy has been 
shown to be crucial in the generation of shape. Indeed 
" a key aspect of shape -petal asymmetry- in the petal lobe 
of Antirrhinum depends on the direction of growth rather 
than regional differences in growth rate" [l9|. To inves- 
tigate the effects of anisotropy, for which our formalism 
is well suited, we study the problem of a free elastic disk 
subject to homogeneous anisotropic growth. 

The model. Since biological soft tissues have a high 
volume fraction of water they are elastically incompress- 
ible; in our notation det A = 1. Furthermore, we assume 
isotropy of the material for simplicity and we define a 
strain energy function W = J2Ts=o ^rs{^i — 3)''(l2 — 3)" 
where Ii and X2 are the principal invariants of the tensor 
A* A. Note that any of the common constitutive relation- 
ships can be described by a series of this form . After 
the deformation process, the sheet, of lateral size L and 
thickness H, is described by the displacement field: u — 
X — X and we define ({X, Y) = uz{X, Y, 0) the displace- 
ment of the middle surface. When the growth rates = 
Gij — 6ij are in the range (H/L)'^ < gij <C 1, the scaling 
of the induced strains falls inside the domain of validity 
of the FvK model. Since 7? <C i, we also apply the mem- 
brane assumption that states axz = cfyz = czz = 0. 
Under these assumptions, after a rather long calcula- 
tion, we derive two main conclusions: fzj The in-plane 
Green tensor, defined as E = ^ (A*A — I) can be ad- 
ditively decomposed into geometric and growth-induced 
part: E = i (F*F — G*G). (ii) All materials behave ac- 
cording to the constitutive equation cr = ^ (E — pi), p 
being the hydrostatic pressure associated with the incom- 
pressibility constraint and Y — 6(coi -I- cio) the instan- 
taneous Young modulus. Thus all thin elastic samples, 
undergoing small (but finite) deflexions follow a gener- 
alised Hooke's law whatever the constitutive relationship 
is, as previously noted in f2]|. The pressure is given 



by the assumptions azz = 0, which imphes p = Ezz- 
Moreover, Ezz can be expressed in terms of the in-plane 
components of the Green tensor, since the incompress- 
ibility constraint yields TrE = 0. Once these results are 
established, we derive the equilibrium equations using 
the principle of minimal energy. They can be written in 
terms of the off-plane displacement and stresses: 
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where Einstein summation convention is used, indices run 
from 1 to 2, _D = YH^/9 is the bending rigidity of the 
plate and </> is a source of mean curvature linked to the 
growth tensor via <j) = Div(GG*) • ez- Using the Airy 
potential x- 

D{A^C~A<j>)-2H[x,C] = P, (2a) 
A\ + Ei[C,C]-i^) = 0, (2b) 

where the [.,.] operator is defined in [22] and the func- 
tion ip appea ring in ([2]) is a source of Gaussian curvature 
explicited in (23|. 

The sets of equilibrium equations ([T]) and ([2]) are a gen- 
eralization of the well known FvK theory of thin plates, 
to which they reduce in absence of growth, i.e, G = I. 
First consider the case where there exists a bijection u 
such that G = Divu. Then, provided that u is consis- 
tent with the boundary conditions, then u is a solution of 
the system of equations with zero elastic energy. In this 
case, the growth process is compatible (with itself and the 
boundary conditions) and there is no need for an elastic 
process, the shape being entirely determined by growth. 
Needless to say, this case is atypical and in general growth 
is incompatible, so neither G nor A are gradients of a de- 
formation field (a one-to-one mapping). Thus the grown 
" state" cannot be physically achieved and is not referred 
to as a configuration. For large deformations however 
(C 3> H), the problem can be simplified. Indeed the 
bending term D {A^( — Acf)) can be neglected and a so- 
lution that cancels the in-plane stresses is a solution of 
[Cj C] = V': called a Monge- Ampere equation. Once this 
equation is solved, the parameters appearing in this solu- 
tion can be selected through minimization of the bending 
energy. For moderate deflections, i.e. ( ^ H, both bend- 
ing and stretching terms are of the same order and the 
solution of zero energy is a surface with prescribed cur- 
vatures, which does not always exist; for example there 
are no surfaces that has positive Gaussian curvature and 
zero mean curvature. It is known that inhomogeneous 
growth can lead to sophisticated surface geometries [l^l , 
so we focus the research on anisotropic growth, which has 
been much less studied. 



3 



The free disk. Consider a disc, of initial radius Ri, 
subject to anisotropic homogeneous growth, with free 
boundaries and no external loading. Referring to a cylin- 
drical system of coordinates {R, Q, Z), the growth tensor 
is diagonal and homogeneous: 



1 + .gi 
G = I 1 + <72 
1 

neglecting the thickening of the plate. If gi and 32, re- 
spectively the radial and circumferential components of 
the growth process, are equal, then growth is homoge- 
neous and isotropic and no residual stress appears: the 
disk remains flat. The relevant control parameter is 
k = 52 — 51- The first case to consider is for ^ ^ ^ 
that induces an off-plane displacement C much smaller 
than H and is thus outside the scope of the present the- 
ory. When k is of order which leads to ( ^ H, all the 
contributions are of the same order and a linear stability 
analysis is performed. We look for a solution in which 
the in-plane fields (displacements Ur, Uq and stresses 
f flfi, (J Re and (Tee) are independent of Q. The off-plane 
displacement, however, can depend on Q. Since the disk 
is free, the boundary conditions imply that there is no 
tension or torque at the free edge and reads <JRR{Ri) = 
aRe{Ri) = 0. The only convergent solution of (jlbp . that 
fulfills these boundary conditions, is (Jrr — (trq = 
leading to Ur{R) ^ ^(52/2 + 51), C/e(i?) = and a 
non-zero hoop stress aee = Assuming a solution 

with discrete axial symmetry: (^{p,Q) = ^(p)cos(m0), 
(with p — R/Ri) equation ([Ta| reduces to 



1 -t- 2m^ 



^(4) + 1^(3) 

P P 

-f ap^ (]^^ — Am? 



K2) 
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m ap 



where a = 



&kR' 



e = (3) 



is a control parameter and is the 



ith derivative of ^ with respect to p. At the free edge 
(p = 1) the zero-torque conditions are not affected by 
the growth process and are described in [23|. To avoid 
singularities at p = 0, we impose ^(0) = and ^'(0) = 0. 
These boundary conditions, together with equation ([3]), 
form an eigenvalue problem for the threshold a. Using 
Frobenius' method [241 , we find the four eigenfunctions 
corresponding to each m. The most unstable mode, oc- 
curring when growth is mainly circumferential (a > 0), 
is characterised by m = 2 -a saddle shape- with a thresh- 
old value of a = 3.08. An axially symmetric solution, i.e 
m = 0, appears when radial growth dominates (a < 0), 
at the threshold value a = —7.82. This simple model 
explains surprisingly well the changes of cap shape that 
the algae Acetabularia acetabulum undergoes during its 
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FIG. 2: Top: The two first destabilized modes, (a) On the left 
fc > 0, the disc adopts a saddle shape, with m = 2, at the thresh- 
old a = 3.08. (b) On the right fc < and the disc adopts an 
axially symmetric shape characterized by m = 0, at the threshold 
a = —7.82. Bottom: shape changes in the Acetabularia algae, the 
figures indicate the fraction of algae that undergo the shape tran- 
sition from an initial population of 85 plants, picture drawn from 

m. 



development. Experiments performed in ^25| | show that 
radial growth occurs in the earliest stage of the devel- 
opment, which leads to a symmetric conical shape. At 
later stage however circumferential growth predominates 
to produces the saddle shape (see FIG. 2). 



We now consider large deformations: 1 ^ fc » 



for which C <C H. Since the stretching contribution is 
much bigger than the bending energy, we first solve the 
Monge- Ampere equation [C, C] — V' in which -0 is given 
by -i/; = kS{p)/p in our case. The general solution is a 
cone that has zero Gaussian curvature except at the tip 
of the cone where the effect of bending becomes impor- 
tant and which would require a more precise treatment 
0, 27 1 . We only focus on the outer solution. The equa- 
tion of the cone is simply ({p,<d) — pg{Q). Using this 
expression, the condition that the Airy potential vanishes 
everywhere (so (Pa|) and ((2b)) are satisfied) gives for the 
in-plane displacement field: 



Uo 



5(e)2 Ue 



{g{ef-g'{er + 2k) de 

(4) 



Periodicity in the orthoradial displacement implies 
C/e(PjO) = Uq{p,2-k). Let 5(6) be represented by its 

Fourier series: 5(6) — X^^^o ("n'^^" 
odicity condition leads to: 



1"® + a*e-i"®l Peri- 



00 



{anal{n^ - 1)) = (oq + aS)^ + 2k 



(5) 
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as dessication. The treatment presented in this letter 
makes the inclusion of anisotropic effects and spatial in- 
homogeneities of the growth process easy. We have shown 
that anisotropic growth induces rich structures like curl- 
ing and crumpling. 




FIG. 3: Top: The two modes of minimal energy, far from the tip. 
(a) fc > 0, the shape is a cone with two oscillations (b) fc < the 
cone of revolution is a solution. Bottom: The resultant shapes built 
out of paper {k = 0.25) are in agreement with the prediction. 



A cone of revolution (for which ao is the only non-zero 
coefficient of the Fourier series) can satisfy the periodicity 
condition only if fc < 0; that is when radial growth domi- 
nates. Infinitely many solutions satisfy the condition ^ 
but the selected shape must have minimal bending en- 
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28l |. The bending contribution reads: 



2 5^(a„a*("'-l)') + (ao+a^)' 



(6) 



Calculating the sum in ^ we find the solutions of 
minimal outer bending energy: 



fc<0 C(^,B) 



(7a) 



fc>0 C(-R,©) i?W^cos2e (7b) 

V o 

For large deformations, those predictions can be easily 
checked by constructing a cone from a disc of paper in 
which a sector defined by two radii is withdrawn and then 
either replaced by a bigger one or just glued to close 
it. This simple demonstration illustrates the fact that 
singularities can arise from growth as observed in dead 
leaves or in the leaves infected by the CLCrV. 

Conclusion Using the formalism introduced by Ro- 
driguez et al, we have developed a theory describing the 
behaviour of thin elastic bodies subject to growth. By 
expliciting the sheet's small thickness, we showed all ma- 
terials behave according to a generalized Hooke's law and 
the equilibrium equations generalize the FvK equations 
accounting for growth. This extension describes a broad 
range of physical phenomena involving mass reorganiza- 
tion, from biological growth to thermal dilatation, as well 
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